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Superconducting LaOFeAs family of even electron metallic systems have striking resemblance to 
odd electron cuprates. We suggest that this resemblance is caused by presence of two coupled 2D 
resonating valence bond systems in LaOFeAs . Bond charge repulsion and Hund coupling fuse the 2 
species of valence bonds into closed strings with a Haldane gap, resulting in a quantum string liquid. 
A pair of doped holes (electrons) creates an open string and remain at the ends as holon (doublon). 
Charge ± 2e singlet strings condense to produce high Tc superconductivity. Higher Tc's are likely 
in our string route to High Tc superconductivity, when competing orders are taken care of. 



Progress in material science and technology continues 
to generate novel many body quantum states. In the 
past decades we have witnessed the birth of quantum 
Hall states, heavy fermions, RVB superconductors, quan- 
tum spin liquids, orbital liquids etc. In the wake of these 
progress, a wealth of novel ideas and theoretical tech- 
niques, with far reaching implication to basic science and 
technology have emerged. Recently superconductivity in 
LaOFeP with a T c as 5 K was discovered PQ, followed by 
an exciting discovery [2] of superconductivity with a T c w 
28 K in electron doped LaOi_ x F x FeAs . Following heels, 
hole doped Lai_ x Sr x OFeAs has been synthesised [5] with 
a T c w 25 K, making electron and hole doping nearly 
symmetric, reminding us of particle-hole symmetric one 
band Hubbard model physics. More recently |H re- 
placement of La by other rare earths Sm, Ce, Nd and Pr 
have resulted in a substantial increase in Tc, 52 K being 
the highest so far. 

The even electron character of Fe 2+ in LaOFeAs , 
orbital degeneracy and observed large superconducting 
T c makes it an unusual high Tc superconductor. Is na- 
ture unfolding another novel organisation of strongly cor- 
related electrons in this 3d 6 system ? In this letter we 
suggest formation of a novel quantum string liquid (QSL). 
The beads of the strings are two types of valence bonds of 
two spin-i resonating valence bond (RVB) systems avail- 
able at the fermi level in LaOFeAs . Fusion of valence 
bonds into strings arise from mutual bond charge repul- 
sion of two different types of overlapping valence bonds 
and Hund coupling. The closed singlet strings have a 
Haldane type correlation and a spin gap [6 

After a brief introduction we write down a model and 
present an RVB mean field theory. Going beyond mean- 
field theory involving 3 different projections brings out 
the string content of the new superconductors. Our es- 
timate of superconducting Tc gives a large value in the 
range 100 to 200 K. We then briefly indicate how the 
underlying string structure can bring about real space 
orders that will in general compete with superconductiv- 
ity. 

The layered transition metal monopnictides , LaOTX 
(T = Mn, Fe, Co Ni; and X = P, As) exhibit a variety of 



low temperature phases: LaOFeP and LaONiP are both 
superconductors |7 with a low T c w 4 and 6 K. LaOCoP is 
a ferromagnetic metal with a T c ~ 50 K and LaOMnP is 
an insulator [2J. Further, iron monopnictides offer some 
surprise: stoichiometric LaOFeP is a superconductor, 
whereas LaOFeAs becomes a superconductor only after 
doping. The spin susceptibility of LaOFeAs is Pauli like 
but considerably enhanced compared to cuprates. There 
is a resistivity anomaly in LaOFeAs at about 150 K. 
Transport, specific heat, H c2 and tunnelling measure- 
ments have been reported. The hall coefficients of elec- 
tron doped LaOi_ x F x FeAs is negative[S]. In an exciting 
work [3], Sr doped Lai_ x Sr x OFeAs has been synthesised, 
resulting in a hole doped superconductor with a T c w 25 
K. The hall constant of Lai_ x Sr x OFeAs is positive. The 
superconducting Tc has steadily increased from 26 K all 
the way up to « 52 K in doped PrOFeAs and NdOFeAs . 
Recent experiments have also confirmed a suggested spin 
density wave order [9 . 

From theory point of view, the first electronic struc- 
ture calculation was performed by Lebegue[T0] for su- 
perconducting LaOFeP . Following the discovery of 
high T c superconductivity in LaOi_ x F x FeAs , several 
groups have performed LDA and more sophisticated 
calculations |1 lj for undoped and doped LaOFeAs and 
also have discussed the importance of strong electron 
interactions in this d-band system. All these calcula- 
tions suggest a dominantly non-phononic mechanism of 
superconductivity. Nature of Spin density wave order in 
LaOFeAs has also been discussed. Different mechanisms 
and order parameters have been also proposed for the 
superconducting state |12j. 

LaOFeAs is a layered material with alternate stack- 
ing of LaO and FeAs layers. For modelling purposes it 
is important to remember that Fe atoms form a square 
lattice, a quasi 2 dimensional metal. Each Fe atoms is 
tctrahedrally coordinated by As atoms. 

Fe 2+ in LaOFeAs is in 3d 6 configuration. LDA calcula- 
tions show 5 overlapping 3d bands spread in the energy 
range -2 to +2 eV. One sees three filled bands and 2 
empty bands, that have a small overlap. The overlap 
leads to small fermi surface pockets with fermi energy 
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~ 0.2 eV. We suggest that this is not a good starting 
point to build low energy physics; we have to take care 
of strong correlation effect. We find that strong corre- 
lations splits the 6 electron LaOFeAs bands into a pair 
of spin- 1 resonating valence bond systems close to the 
fermi level, as shown schematically in figure 1. We call 
it as 2 RVB system. We will present this analysis in a 
separate communication. 
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FIG. 1: Schematic picture of splitting off of two half filled 
bands at the fermi energy, from filled bands as coulomb re- 
pulsion is increased. 

Our starting point of this paper in discussing 
mechanism and theory of superconductivity in doped 
LaOFeAs will be a 2 RVB system. A major prediction of 
our model building is that ARPES should see two large 
fermi surfaces close to half filling, rather than small fermi 
surface pockets predicted by LDA calculations. 

LDA calculations finds that the tetrahedral crystal 
field splitting and a large direct 3d orbital overlap spread 
the spectral weight of every 3d orbital over the entire 
band width of 4 eV. As for the symmetry of the Wannier 
orbitals of these two bands are concerned, band order 
suggested by tetrahedral crystal field splitting need not 
be a good guidance, because direct Fe 3d orbital overlap 
produce orbital mixing and compete with crystal field 
splitting. Based on symmetry grounds Wannier orbitals 
will have the same symmetry property of the 3d orbitals; 
however, they will have strong admixture with various d 
orbitals of neighboring sites. 

Assuming that the metallic state maintains the square 
lattice symmetry, there are different possibilities for 
choosing symmetry of the two Wannier orbitals: i) a pair 
from 3d x 2_ y 2 , 3d xy , 3d z 2 , whose \ip\ have square pla- 
nar symmetry' or ii) 3d zx and 3d zy . The second choice 
will give 2 spatially anisotropic bands, but together will 
maintain the square lattice symmetry. 

So we choose a fairly general two orbital Hamiltonian. 
Key parameters of our model arc inter and intra orbital 
Hubbard U, ~ 3 to 4 eV, bond charge repulsion V12 ~ 
1 eV and Hund coupling ~ 0.7 - 1.0 eV. Width of an 
individual 3d band is about 2 eV, giving a hopping pa- 
rameter between nearest neighbour Wannier orbitals, t 



~ 0.5 eV. 

The Hamiltonian of our 2 RVB systems is a two orbital 
Hubbard model: 

+ V 12 n ijl^j2 ~ J H ^a^WlC^) 2 (1) 

(ij) i M 

Here /x, v = 1,2 represents Wannier orbitals. Hopping 
is assumed to exist only among same type of nearest 
neighbour orbitals. Second line of the above equation 
couples the two systems, through bond charge repulsion 
and Hund coupling. The operator n^ = ^J2a( c lfj,cr + 

c jfj,a)( c incr + c jn<j) counts number of electrons in bond- 
ing state of ^t-th orbital connecting neighboring sites 
i and j. The bond charge repulsion is a coulomb in- 
teraction term,Vi2 = / \ipa{r) + ipj^r)] 2 j^pj\tp l2 (r') + 
ipj2(r')\ 2 drdr' , where are the two Wannier orbitals 
at site i. 

Starting from the above Hamiltonian, using a superex- 
change perturbation theory, we derive the following ef- 
fective Hamiltonian for the optimally doped case. It is a 
sum of two t-J models: 



H cS = H U1 + H U2 = -^2 hjAn^j^ + h - c - + 
- ^ ' Sj M — -n ltl n jtJi ) (2) 



with three local constraints: i) + 7^ or 2, for 
electron or hole doped cases, ii) bjj^&y^ 7^ 2 and iii) 

(E P i^^ c ¥ / ^ °- The superexchange J~ j. At 
the Hamiltonian level the 2 RVB systems are decoupled; 
however the three constraints couple two RVB systems in 
a non trivial fashion. The first, double/zero occupancy 
constraint is well known in t-J model. The second and 
third are new in the context of superconductivity theory. 
The second one tells us that two neighboring sites con- 
taining 2 electrons each cannot form two covalent bonds, 
because of the bond charge repulsion. The third con- 
straint is dictated by Hund coupling, which favours max- 
imal spin at a given site. 

We will discuss the above model Hamiltonian, using 
RVB theory approach [13j [14]: i) solve the unconstrained 
Hamiltonian in a mean field theory and ii) perform all 
projections in the resulting wave function, for further 
analysis. The unconstrained Hamiltonian is the same as 
the one solved in the first RVB mean field theory [14] of 
cuprates. Later several important improvements [i~5] have 
been made on that approach. We can use all those re- 
sults. Formally the wave function we wish to analyse for 
superconductivity in our 2 RVB system (for hole density 
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x) is the following: 
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erator. Further, Pq is the usual Gutzwiller projection 
which prevents double occupancy in any of the orbitals. 
P b avoids avoids two singlet bonds connecting two neigh- 
boring sites. P$ projects out the singlet spin component 
of two electrons at any site. The pair function <pij is the 
variational wave function. For a single RVB system, in a 
square lattice, the above function has d x 2_ y 2 symmetry 
and leads to a superconducting state with nodal quasi 
particles. In the presence of coupling between two RVB 
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FIG. 2: Fusion of Valence Bond States. 2a) VB states of the 
two systems. 2 b) String like organization after fusion. 

systems key questions are: i) do the three projections pre- 
serve d-wave superconductivity of a single RVB system? 
and ii) are there going to be some new physics? These are 
difficult questions to answer. However, the available phe- 
nomenology suggests superconductivity, similar to that 
suggested by 2D repulsive Hubbard model close to half 
filling survives. In what follows, we will address these 
question from the point of view of short range RVB states 
and find that projections bring an entirely new quantum 
liquid namely quantum string liquid and associated rich 
possibilities. 

We also develop a quantum mechanical basis of sin- 
glet string states, by an appropriate fusion of two short 
range RVB states. We will start with the case of no dop- 
ing. Consider two nearest neighbour valence bond (VB) 
states shown in figure 2a. The VB's of two systems are 
denoted by dark and shaded bonds. If we choose two 
arbitrary VB states and fuse them, we will get overlap- 
ping bonds, in general. Such states are energetically not 



favourable because of bond charge repulsion energy V12. 
The two valence bond states in figure 2a do not have 
bond overlap. When we fuse them we get a state shown 
in 2b. From simple topological considerations it follows 
that the resulting state is a set of closed strings. Open 
strings, if they appear, carry spinons at their ends. Also 
every configuration of closed and open strings (made of 
nearest neighbour bonds) that fill the lattice can be ob- 
tained by a fusion of two unique VB states. 

Short range VB states form a overcomplete set of states 
to describe the physics of a spin-^ Heisenberg model. The 
overcompleteness make different VB states linearly inde- 
pendent and not orthogonal. As our fused states are di- 
rect product of two RVB states, we can use known results 
of single VB overlap properties to study our combined 
system. It should be remembered that the total number 
of string states is simply not square of the number of 
possible valence bond states, because of the constraints. 

Let us calculate energy expectation value of the bond 
repulsion and Hund coupling term in a string state. By 
construction bond repulsion energy is zero. Since we have 
two valence bonds meeting at every site, the total spin 
value at a given site is fluctuating: with probability | 
it has value zero and with probability | it has value 1. 
Thus the average value of the Hund coupling energy is 
|th of the maximum possible value 2J#. It is possible to 
gain the extra Hund coupling energy |2J# projecting out 
the singlet component at every site. Or one can replace 
every string state formally by the exact ground state of 
the Haldane gapped nearest neighbor Heisenberg anti- 
ferromagnetic spin-1 chains. Our strings already have an 
orbital order as the beads alternate along the chain. Con- 
verting them into Haldane chains brings an additional 
topological order [IS], 



charge j-2e 
cooper/pair 



FIG. 3: charge -2e Cooper pair, a spin-1 spinon pair and an 
electron, as different string states. 

Now we briefly discuss topological excitations in our 2 
RVB system. Picture 2b already shows how a two spinon 
(spin-1) state appears as an open string. As every un- 
paired spin or doublon or holon can occur in one of the 
two orbitals of a given site, they carry an orbital quantm 
number as well. A bound spinon pair (figure 3) is likely to 
be a low energy excitation in our metallic system. When 
we add an electron to the insulating valence bond state, 
we get a spinon-doublon composite as shown in figure 3. 
When two electrons are added, it is energetically advan- 
tageous to get rid of two unpaired spins. That is, the 
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unpaired spin will disappear as singlets and we will be 
left with one open string with two doublons (charge -e) 
at the ends (figure 3). 

Using our fused short range VB states or equivalently 
the string state basis, we can describe superconducting 
ground states and excited states. What is interesting is 
that the fusion with the constraints has resulted in an 
incompressible liquid of strings. There are interesting 
questions about how strings quantum fluctuate in space 
and manage to create a string liquid. This is a non triv- 
ial question, as the string states are also non-orthogonal. 
Even the simple question of evaluation of overlap of 2 
string states becomes complicated, if we keep the strings 
in Haldane gapped state. More remains to be done to 
take advantage of the string structure of our 2 coupled 
RVB states. The nature of holon as well as spinons at 
the end of the open strings are very interesting. It is 
known that open Haldane chains support spin-| excita- 
tions at the boundary 1 17], However, the gap moment is 
not strictly localised on the chain end. The excess spin 
density oscillates and decays exponentially, giving a char- 
acteristic size for the spinon. We expect the same for the 
holon as well. 

An open string with charges holon/doublon at both 
the ends is our cooper pair. It is not obvious how it 
will modify the nature of superconducting state. While 
the overall ODLRO and phenomenology may resemble 
the standard BCS superconductor, there may be subtle 
topological orders and non trivial excitations in our QSL 
superconductor. It needs to be explored. 

String structure suggests, depending on the size of the 
open string (dictated by factors such as resonance en- 
ergy, coulomb repulsion etc.) a binding mechanism for 
our holon pairs. Is this an additional pairing energy ? 
The final superconducting state will be a coherent super- 
position of the resonating charge string configurations. A 
key parameter that determines the superconducting Tc 
will be the effective mass of the 2e open string. That will 
give us a Kosterlitz-Thouless type of scale ksT c w jrrT' 
where n is the carrier density per unit area. This ex- 
pression is very similar to the expression or Tc in RVB 
theory, suggested by the condenseation of charge valence 
bonds or holons. Thus we expect a maximum Tc in the 
range of 160 to 200 K, perhaps exceeding cuprates. 

Since we have string like entities, there will be a ten- 
dency for them to have liquid crystalline type order, spin 
order and charge order, encouraged by unscreened long 
range interaction at low doping and electron lattice cou- 
pling. Such real space organization will in general reduce 
superconducting Tc. These are competing phases, very 
much like in cuprates. If one can engineer materials, such 
as Tl or Hg multilayer cuprates, where charge and spin 
order tendencies are supressed, superconducting Tc's can 
go higher than 52 K that has been observed so far. 

Experimentally one sees a metallic state in undoped 
LaOFeAs . This is likely to arise from two possibilities: 



i) the two reference spin- ^ Mott insulators having an oc- 
cupancy of 1 + xo and 1- xo per site; i.e., they are self 
doped (xo << 1) for energetic reasons and ii) one refer- 
ence state remains insulating whereas the other remains 
strongly correlated but metallic. In the second case a 
commensurate long range order at (tt, it) can easily oc- 
cur in the metallic band. External doping may selective 
enter the Mott insulator and we get superconductivity, 
somewhat contaminated by the first metallic band. 

The spin susceptibility of the undoped metal 
LaOFeAs and LaOi_ x F x FeAs are both high compared 
to the metallic cuprates. This is likely to arise from the 
presence of low energy spin-1 excitations in the quantum 
fluctuating closed strings. 

It is a pleasure to acknowledge Venky Venkatesan, who 
brought to my attention reference 2, and knowing my 
fascination for new superconductors, called it a gift. 

Postscript: As we were completing this manuscript a 
paper by Tao Li (cond-mat/0804.0536) appeared. It in- 
troduces a model of two strongly correlated 3d zx and 
3d zy bands and presents a mechanism of spin density 
wave and superconductivity. 
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